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Stellar Relaxation Processes Near the Galactic 
Massive Black Holef 

By Tal Alexander^ 2 

^Faculty of Physics, Weizmann Institute of Science, PO box 26, Rehovot 76100, Israel 
■^William Z. and Eda Bess Novick Career Development Chair 

The massive black hole (MBH) in the Galactic Center and the stars around it form a unique stellar dynamics 
laboratory for studying how relaxation processes affect the distribution of stars and compact remnants and 
lead to close interactions between them and the MBH. Recent theoretical studies suggest that processes 
beyond "minimal" two-body relaxation may operate and even dominate relaxation and its consequences 
in the Galactic Center. I describe loss-cone refilling by massive perturbers, strong mass segregation and 
resonant relaxation; review observational evidence that these processes play a role in the Galactic Center; 
and discuss some cosmic implications for the rates of gravitational wave emission events from compact 
remnants inspiraling into MBHs, and the coalescence timescales of binary MBHs. 



1. Introduction 

The M.~ 4x10*^ massive black hole (MBH) in the Galactic Center (GC) dEisenhauer et al. 



20051: iGhez et a/.ll2005h and the stars around it are the closest and observationally most acces- 
sible of such systems. Observations of the GC thus offer a unique opportunity to study in great 
detail the effects of the MBH and its extreme environment on star formation, stellar evolution 
and stellar dynamics, and the interactions of stars and compact remnants with the MBH. 

Here the focus is stellar relaxation processes. Relaxation plays an important role in a wide 
range of phenomena that involve close interactions with a MBH (the "loss-cone problem", ^l.ll i. 
Such processes include gravitational wave (GW) emission by compa ct remnants inspiraling into 
a MBH ("extreme-mass ratio inspiral ev ents" (EMRIs, see rev iew bv lAmaro-Seoane et al. ]l2007l 
tidal flar es from tidal disruption events ([Fra nk & Rees' '1976'), tidal capture and tidal scattering 
of stars ([Alexander & Morrisll2003t [Alexander & Liviai200L) . 3-body exchanges with binaries 
leading to the capture of stars on ti ght orbits ar ound the MBH and the ejection of hyper veloc- 
ity stars (HVSs) out of the galaxy (lHillslll988b . the orbital decay and coalescence of a b inary 
MBHs (and the "last parsec stalling problem", see review by Merritt & Milosavhevicl2005 ), and 
perhaps also the formation of ultra-luminous X-ray sources (UL Xs) in star clusters f ollowing 
stellar capture around an intermediate mass black hole (IMBH) (iHopman et a/.ll2004l) . Relax- 
ation processes are possibly linked to the presence and properties of unusual stellar populations 
that are observed near MBHs, such as the c entral "S-star" cluster, the s tellar disks in the G C 
dEisenhauer et al. I2OO5I: iPaumard et al 1120061) and the stellar disk in M3 1 ( iBender et al. ]l2005h . 

Dynamical relaxation by star-star interactions is inherent to the discreteness of stellar systems. 
In the absence of additional mechanisms to randomize stars in phase-space, standard 2-body 
stellar relaxation assures a minimal degree of randomization, albeit one that could be too slow 
to be of practical interest. This review will discuss relaxation processes beyond standard stellar 
relaxation, which operate on much shorter timescales, or else operate in a qualitatively different 
way: massive perturbers (©, strong mass segregation (© and resonant relaxation (® . 

The dynamical properties of the GC, specifically its short 2-body relaxation time and high 
stellar density, are probably not typical of galaxies in general 01.2I) . However, dynamical pro- 
cesses that can be probed by GC observations have implications beyond the GC. In particular, the 
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Milky Way is the archetype of the subset of galaxies with low-mass MBHs that are key targets 
for planned space-borne gravitational wave detectors, such as the Laser Interferometer Space 
Antenna (LISA). GC studies may help understand the effect of such relaxation processes on the 
open questions of the cosmic EMRJ event rate and the EMRI orbital characteristics. 

Before turning to a discussion of the non-standard relaxation processes that are expected to 
operate in the GC, it is useful to briefly review the dynamics leading to close interactions with a 
MBH (loss-cone theory) and the dynamical conditions in the GC. 

Infall and inspiral into a MBH 

Stars can fall into the MBH either by losing orbital energy, so that the orbit shrinks down to 
the size of the last stable circular orbit (r-Lsco = 3rs for a non-rotating MBH, where the event 
horizon is at the Schwarzschild radius = 2GM,/c^ ), or by losing orbital angular momentum 
so that the orbit becomes nearly radial and unstable (periapse rp < 2rs for a star with zero orbital 
energy falling into a non-rotating MBHjJ. The timescale to lose energy by 2-body scattering, 
Te = \E/E\ is of the order of the relaxation time, 

Te^Tr^ {M,/M,f Tdyn{r)/N4< r) logN,{< r) , (LI) 

where A'^^(< r) is the number of stars inside r, Tdy„(r) ^ ^Jr^/GM, is the dynamical time 
and spherical symmetry and a Keplerian velocity dispersion are assumed, ^ GM,/r. The 
maximal angular momentum available for an orbit with energy E is that of a circular orbit, 
Jc{E) = GMtj \J2E (using here the stellar dynamical sign convention E = — — 0(r) > 0). 
The timescale for losing angular momentum, Tj = \ J I J\, can be much shorter than Te when 
J < Jc, since 



Tj = \JIJ,(E)fTE. (1.2) 

As a consequence, almost all stars that reach the MBH, and are ultimately destroyed by a close 
interaction with it, do so by being scattered to low- J "loss-cone" orbits (near radial orbits with 
J < Jic ':^ \/2GM»Q, where q is the maximal periapse requ ired for the close interaction of interest 
to occur Iprank ARees"l976; 'Lightm an & Shapirdll977l ). The rate of close interaction events, 
Fic, is set by the replenishment rate of stars into the loss-cone. When the replenishment mecha- 
nism is diffusion in phase space by 2-body scattering, Fic oc , which is typically a very low 
rate. Close to the MBH, at high-i?, where the relative size of the loss-cone in phase-space is large 
( J/c/ Jc \fqE), relaxation is too slow to replenish the lost stars, and the loss-cone is on average 
empty. Farther out, at low E, where the loss-cone is small, relaxation can replenish the lost stars, 
the loss-cone is full (isotropic distribution of stars) and the local replenishment rate is maximal. 
Nevertheless, the contribution to the total replenishment rate from the low-i?, full loss-cone re- 
gions of phase-space, where the timescales are longer and the stel lar densities lower, remain s 
small compared to that from the empty loss-cone regions at high-£^ (ILightman & Shapiroll977t ). 



The observational and theoretical interest in such close interactions motivated numerous inves- 
tigations of alternative efficient loss-cone replenishment mechanism, suc h as 2-bod y relaxation in 



non-spherically symmetric potentials (.Magorrian & Tremaine 1999: Berczik et a/.l2006.). chaotic 

004 HGerhard & Binney| l985[ 
perturbers fehao et al.' '2002' 



age 

orbits in triaxial potentials dNorman & Silkll983l:lMerritt & Poonl2004l:lGerhard & Binney|l9"85| 



massive 



Hollev-Bockelmann & Sigurdsson' '2006'), relaxation by 
Perets ef al. 2007). resonant relaxation (Ra uch & Trema ine 1996; Ra uch & Ing alls 1998: Hopman & Alexandei 
2006d:lLevinl2007h. or perturbations by a massive accretion disk or a secondary IMBH (IPolnarev & ReesI 



1994l:lLevin et al I2OO5L 



f If the stars are tidally disrupted before falling in the MBH, the relevant distance scale is the tidal 
disruption radius rt ~ ii*(M. /M*)^'^'* > rs rather than the event horizon ra. 
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Close interactions with a MBH fall in two dynamical categories ([Alexander & Hopmanl200i : 



infall processes, such as tidal disruption, where the star is destroyed promptly on its first close 
encounter with the MBH, and inspiral processes, such as GW EMRI events, where multiple 
consecutive close encounters are required for the orbit to gradually decay. The infall takes about 
an orbital period, the time to fall from the point of deflection to the center, whereas the inspiral 
process takes much longer, depending on the energy extraction efficiency of the dissipational 
mechanism involved (for example GW emission, tidal heating or drag against a massive accretion 
disk). In most cases the dissipated energy is a steeply decreasing function of the the periapseQ 
and so the inspiral time scales with the number of periapse passages, and hence with the initial 
orbital period. 

An infall or inspiral event can occur only if the star, once deflected into the loss-cone, avoids 
being re-scattered out of it (and in the case of inspiral, also avoids being scattered directly into the 
MBH). Because inspiral processes are slow, stars can avoid re-scattering, complete the inspiral 
and decay to an interesting, very short period orbit with high emitted dissipative power, only if 
they are deflected into the loss cone from an initially short period orbit, with E > EcHt- Figure (HJ 
shows a schematic description of the phase-space evolution of infalling and inspiraling stars, and 
the emergence of a critical energy scale. For inspiral by GW emission into a A/, ^0(10^ ^^o) 
MBH, E'ciit corresponds to an initial distance scale of rciit ^ 0.01 pc (the ansatz r ^ E = 
GM,/2a, is assume d here, where a is the Keplerian semi-major axis). The EMRI event rate is 
then approximately dHopman & Alexand er 2005) 



r 



' iVGw(< ?'crit)/Jfl(rcrit) OC A^Gw(< 7'crit)A^*(< '"crit )/Tdyn (T-crit ) , (1-3) 

where AfGw(< f) is the number of potential GW sources (compact remnants) within distance 
r of the MBH. A critical energy can be similarly defined for infall processes. Because infall is 
much faster than inspiral, i?ciit is much lower ( rerit much larger). For example, the critical radius 



for tidal disruption in the GC is rcrit ~ fcwpc ( Lightman & Shapiro|[l977 : Is ver&Ulme3l 19991: 



Magorrian & Tremainelll999l) . Most of the stars that infall or inspiral originate near 



^crit ■ 



Equation ( 11.31 ) shows that the degree of central concentration of compact remnants strongly 
affects the EMRI event rate. Mass segregation therefore substantially increases the predicted 
EMRI event rate from inspiraling 0(10 MfT^) stellar black holes (SBHs), w hich are the most 
massive, long-lived objects in the population jHopman & Alexandei 2006^; ®. Similarly, the 



capture of compact remnants very near the MBH by 3-body exchanges between the MBH and 
binaries (© can also strongly affect the EMRI rate (Perets, Hopman & Alexander, 2007, in 
prep.). The dependence of Tie on Tn is not trivial, since rent itself depends on Tf>: the shorter 
the relaxation time, the faster stars are scattered into the loss-cone, but also out of it. Detailed 
analysis shows that the two effects cancel out for ocr^'^/^ stellar cusps. Since in most galactic 
nuclei the logarithmic slope of the density profile is not much di fferent from —3/2, the EMR I 
rate is expected to be roughly independent of the relaxation time dHopman & Alexandej|2005l) . 



It should be emphasized that this result applies only to 2-body relaxation, and needs to be re- 
examined if other loss-cone replenishment mechanisms dominate the dynamics. 

1.2. The dynamical state of the stellar system around the Galactic MBH 

The stellar system around the Galactic MBH is expected to be in a state of dynamical relaxation 
in a high density cusp. This is a direct consequence of the low mass of the Galactic MBH and 
of the M,/cr relation, the tight observed correlation between the mass of central MBHs and 
th e typical velocity dispersion in the bulges of the ir host galaxies, M, oc a^, where 4 < /? < 



5 (iFerrarese & Merritl I2OOOI: iGebhardt et a/.ll200ol) . /3 = 4 is assumed here for simplicity; the 



conclusions below are reinforced if /3 > 4. 



i E.g. the GW energy emitted per orbit scales as AE ix{A'hc^ /M,)(rp/rs) (lPetersll964h . 



4 T. Alexander: Relaxation Processes Near the Galactic MBH 




background" 



logJic log J logJic log J 

Figure 1 . A schematic representation of the phase-space (log E, log J) trajectories leading a star to the 
MBH. Each segment of the random-walk trajectory represents the change in the phase coordinates over 
some fixed time step At. The shaded areas on top {E > -Blsco) and on the left (J < Jic) are regions of phase 
space where no stable orbits exist. The diagonal boundary on the right is s the maximal angular momentum 
Jc{E). Left: Infall without dissipation. A star with initially high J is scattered with roughly equal relative 
magnitude in E and J. Eventually a random kick will send it to a low- J orbit, where J-scattering is much 
faster than i5-scattering, making it plunge directly into the MBH. Right: Inspiral with dissipation. Energy 
dissipation by the emission of GW can lead to very rapid orbital decay on low- J orbits, faster than the mean 
time between scattering events, thus enabling the star to reach a short-period orbit with detectable GW 
emission (narrow horizontal shaded strip on top). Statistically, nearly all stars with initial energy E > -Ecrit 
will ultimately inspiral into the MBH, while nearly all stars with E < Ecrit will ultimately plunge into the 
MBH, following a trajectory similar to the one depicted in the left panel. 



The MBH radius of dynamical influence is conventionally defined as ~ GA/./cr^ oc All^^. 
The mass in stars within the radius of influence is of the order of the mass of the MBH, so their 
number is Nh ^ M,/Mi,, where is the mean stellar mass, and the average stellar density 
within r/i is Uh ~ Nh/rf^- The two-body relaxation time at Vh is Tr ~ (A/,/i\/^) Th/Nh. It then 
foflows that Tr (X mI'^ and fih cx M,^'^. Evaluated for the Galactic MBH, Tr - 0(1 Gyr) < 
tn (the Hubble time) and Uh ^ 0(10^ pc"'^). The short relaxation time implies that the system 
will return to its relaxed steady state following a major perturbation, such as a merger with a 
second MBH (iMerritt & Szellll2006l: iMerritt ef a/.1l2007h . Note that Tr > t,, for a MBH only a 
few times more massive than the Galactic MBH. The GC is thus a member of a relatively smaU 
subset of galaxies with high-density relaxed stellar cusps. 

A relaxed stellar system is expected to settle into a power-law cusp distribution, rii, cx r~", 
(® . The high stellar density in a steeply rising cusp allows star-star and star-MBH interactions 
to become frequent enough to be dynamically relevan t and observationally i nteresting (Eg. II. 3b . 
For example, the rates of both tidal disruption events ( Wang & Merrittll2004 ) and EMRI inspiral 
events ( Hopman & Alexandeill2005 ) scale inversely with the MBH mass, F cx A^/j/Tr cx A/, 



2. Massive perturbers 

2.1. Massive perturbers and the loss-cone 

The relaxation time (Eg. II. Il l is proportional to {AI'^n^,)~^ . This can be readily understood by 
considering the "F ~ nvY," collision rate between stars of mass M^, and mean space density 
in volume V, n^, = Ni,/V, where the cross-section S ^ Trr'l is evaluated for collisions at the 
capture radius Tc ~ IGM-^jv^ , the minimal radius for a soft encounter with a typical velocity v. 
The rate of scattering by stars is then F^, ~ A/^ /v^ ~ ^ (integration over all colHsion radii 
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increases the rate only by a logarithmic Coulomb factor). When the system also contains a few 
very massive objects such as giant molecular clouds (GMCs), stellar clusters, or IMBHs (if such 
exist), these massive perturbers (MPs) of mass Mp ^ Mi, and space density Up ~ Np/V <C rii, 
will scatter stars at the capture radius Tc = G{Mi,+Mp)/v'^ at a rate of T p ^ np{Mi,+ MpY / ~ 
UpMp I v^. MPs could well dominate the relaxation even if they are very rare, as long as 



H2 = M^Np/M^N^>l 



(2.1) 



Efficient relaxation by MPs was first suggested by ISpitzer & Schwarzschild ( 1951 , 1953 ) to 
explain stellar velocities in the Galactic disk. Its relevance for replenishing the loss-cone was 
subsequently investigate d in the context of Sola r system dynamics for the scattering of Oort 
cloud comets to the Sun jHiUslll 98 iLlBailevll 1983b . and more recently as a mechanism for estab- 
lishing the M,/a correlation by fast accretion of stars and dark matter dZhao et a/.ll2002h . Here 
the focus is on the consequences of MPs fo r the replenishmen t of the loss-cone, and the implica- 
tions for stellar population s in the Galaxy dPerets et al.luQQlh . the coalescence of binary MBHs 
( Perets & Alexandeil2007 ) and for the cosmic rates of EMRIs (Perets, Hopman & Alexander, in 
prep.). 

Loss-cone replenishment by MPs can b e described by the standard loss-cone formalism (e.g. 
Young|[l977 ) with only few modifications ( Perets et a/.ll2007 ). The large size of the MPs is taken 
into account by decreasing the Coulomb logarithm accordingly; the orbital averaging of phase- 
space diffusion due to scattering by stars is done incoherently (sum of squares), while for rare 
MPs, where there may be on average less than one scattering events per orbital period, the aver- 
aging is done coherently (square of sums). 

The relative contributions of relaxation by stars and relaxation by MPs to the total loss-cone 
replenishment rate depend on the size of rcrit relative to the spatial distribution of the MPs 
(rciit increases with the loss-cone size, and in the case of inspiral also with the efficiency of the 
dissipative process). MPs are extended objects, which cannot survive in the strong tidal field of 
the MBH (IMBHs could be the one exception). Generally, MPs in galactic centers could also 
be affected by an intense central radiation field, whether the AGN's or the stars', or by outflows 
associated with accretion on a MBH. These processes introduce an inner cutoff 7-mp to the MP 
distribution. A plausible estimate is tmp This is the case in the GC, where the clumpy 

circumnuclear gas ring lies outside the central 1.5 pc, on a scale comparable to r^. The event 
rates of processes such as tidal disruption of single stars (rcrit ~ f /i) or GW EMRI (rcrit <C Vh), 
where stellar relaxation by itself efficiently fills the loss-cone at rcrit < r < vmp, will not be much 
enhanced by additional relaxation due to MPs (the stellar distribution function (DF) cannot be 
more random than isotropic). In contrast, the event rates of processes whose loss-cone is large, 
and which would have remained empty beyond rup in the absence of MPs, can be increased by 
orders of magnitude by the presence of MPs. Most of the enhancement is due to MPs near rup 
dPerets et a/ ...2007). 

Here we consider two processes with large loss-cones, where MPs play an important role: the 
tidal disruption of stellar binaries of total mass A/12 and semi-major axis ai2 that interact with 
the MBH at a distance rp<rt^ ai2(A/,/M i2)^^^, lead ing to the capture of one star around the 
MBH and the ejection of the other as a HVS dHillsll 19881) . and the orbital decay of a binary MBH 



of total mass A/12, mass ratio M2/M1 = Q <1 and semi-major axis ai2 by interactions with stars 
at a distance rp < 0(012) (the "slingshot effect") ( Begelman ef a/. 198Q|). 

2.2. Massive perturbers in the Galactic Center 

MPs in the GC include GMCs, stellar clusters and possibly IMBHs, if these exist. Direct ob- 

5-100 pc scale are 
and a typical size of 



servational evidence (Fig. |2]i indicates that the dominant MPs on the r 
0(100) GMCs in the mass range lO^^-lO^ A/©, with rms mass of - 10^ A/ 



i?p ~ 5 pc (the quoted range includes an order of magnitude uncertainty in the mass determina- 
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Figure 2. The observed MPs in the GC and their effect on the relaxati on time. Left: Th e observed mass 
function of molecular cloud massive perturbers in the GC (adapted from lPerets et al\ i2007b with permission 
from the Astrophysical Journal). Lower (o) and upper (virial) (□) mass estimate s for th e molecular clumps 
in the inner ~ 5 pc, based on the molecular line observations of lChristopher et al] ( l2005h . and lower (A) and 
upper (virial) (v) mass estima t es for the GMCs in the inner ~ 100 pc of the Galaxy, based on the molecular 
line observations of lOka et al.\ ( 120011) . Right: The relaxation time as function of distance from the Galactic 
MBH due to stars alone, the upper (GMCl) and lower (GMC2) mass estimates of the molecular clumps 
and GMCs and due to upper (Clusters 1) and lower (Clusters2) estimates on the number and masses of 
stellar clusters. The sharp transitions at r = 1.5 and 5 pc are artifacts of the non-continuous MP distribution 
assumed here. GMCs dominate the relaxation in the GC. 



tion), jOka et a/.ir200ll:lGusten & Philippll2004 . and on the r-- 1.5-5 pc scale, 0(10) molecular 
clumpjj with ma sses in the range lO'^-l O^ Mq, with rms mass of ~ 10^ A/q and a typical size 



of ~ 0.25 pc (IChristopher et a/.ll2005h . The ~ 10 observed stellar cl usters ([Figer et al. I ll 9991: 



Borissova et al. 2005h may compete with stellar relaxation ( ^rets & Alexandei l2007l) . Com 



pared to the 2 x 10** 1 Mq stars in the central 100 pc (iFiger ef a/.ll2004 . the GMCs are 
expected to decrease the relaxation time by a factor ^2 ~ 50-5 x lO'^ (Eg. 12. 11 1. Figure Q shows 
a more detailed estimate of the local relaxation time for the various molecular cloud models, 
taking into account, among other considerations, the Coulomb factors. The relaxation time is 
indeed substantially decreased, by factors of 10 — 10^ relative to that by stars alone, depending 
on distance from the center, and on the GMC mass estimates. If IMBHs do exist, then the effects 
of accelerated relaxation will be even stronger than predicted here, and probably extend all the 
way to the center. 

2.3. Galactic and cosmic implications 

With stellar relaxation alone, the empty loss-cone region of MBH-binary interactions is large 
(rt (X 012) and extends out to > 100 pc. However, the MPs that exist in the Galaxy on that 
scale accelerate relaxation, efficiently fill the loss-cone, and thus increase the binary disrup- 
tion rate by several orders of magnit ude, making binary disruptions dynamically and obser- 
vationally relevant (IPerets et a/.ll2007h . Such events, which result in the energetic ejection of 
one star, and the capture of the other on a close orbit around t he MBH, have various pos- 
sible implications. Disruptions of binaries by the Galactic MBH (lHillslll988l : IYu & Tremaine 



f The division of a quasi-continuous medium into individual clouds is somewhat arbitrary, since several 
sub-clumps can be identified as a single cloud, depending on the spatial resolution of the observations and 
the adopted definition of a cloud. For a fixed total MP mass, M = MpNp within a region of size r, the 
relaxation time scales with Np as Tr oc {MpNp)~^ = M~'^Np\ the more massive and less numerous 
the clouds, the shorter Tr. The value of Tr. thus depends on the way clouds are counted. Obviously, the 
statistical treatment of relaxation is valid only for Np S> 1 and Rp <C r. 
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Figure 3. A schematic representation of thie mapping of tire initial binary period distribution to the semi- 
major axis of the tidally captured sta r. Left: The bimodal initial period distribution for old white dwarf / 
main sequence binaries (adapted from lWillems & Kolb 2004il . and its subsequent evolution due to GW co- 
alescence (for the shortest periods) and to slow evaporation by field stars at 100 pc and faster evaporation 
at 10 pc. Right: the resulting semi-major axis distribution of the captured stars, due to scattering by stars, 
which occurs on the 0(10 pc) scale, and due to scattering by MPs, which occurs on the 0(100 pc) scale. 



2003 : Gualandris et al ]|2005 : iBromlev et al were suggested to be the origin of the hyper- 



veloc i ty~B-starjl (y > SOOkms ^), observed tens of kpc away from the GC ( Hirsch et al. 



2005t iBrown et a/.ll2005l: lEdelmann et a/.ll2005l: iBrown et a/.ll2006al) . and the origin of the puz- 



zUng "S-stars" dOould & OuillenI 120031: lOinsburg & Loebll2006l) . a cluster of ~ 10-30 main- 
sequence B-stars (4 Mq < M^, < 15 Mq, main sequence lifespan i^, ~ few x 10''-few x 10^ yr) 



on random tight orbits around the MBH in the central few x 0.01 pc (lEisenhauer et a/.l 12005 



Ghez et al bop s'). Compa ct obje cts captured this way could eventually become zero-eccentricity 



inspiral ( Hopman & Alexandeil 



GW sou rces ([M iller e fa/.l 12005), in c ontrast to high-eccentricity sources typical of single-star 



20051). These two classes of sources are expected to emit very 



different gravitational wave-forms. 

Dynamical arguments and simulations show that on average, ^ 0.75 of MBH-binary encoun- 
ters lead to capture, a nd that the mean s emi-major axis of the captured star is related to that of 
the original binary by (lHillslll988l[l99l 



(a) ~ (A/./Mi2)'/'ai2, (2.2) 

which implies a very high initial eccentricity, 1 — e = rt/ (a) = (Mi2/M,)-'^/'^ ~ 0(0.01). The 
tidal capture process can be viewed as a mapping between the properties of field binaries far from 
the MBH, and the orbital properties of the captured stars: wide binaries result in wide captured 
orbits, and vice versa (Fig. O. The mean velocity of the ejected star at infinity (neglecting the 
galactic potential) is 

(vl) ^ V2GA/ff A/yVai2 • (2.3) 
This translates, for example, to Voo ~ 2000 km s^^ for a 2 x 4 Mq B-star binary with ai2 = 0.2 
AU, well above the escape velocity from the Galaxy. 

Of particular interest is the connection between the HVSs and the S-stars that is implied by the 
binary tidal disruption scenario. The stellar binary m ass ratio distribution is peaked around ^ 1 
( Duquennov & MavoiJl991 : Kobulnickv et al. Il2006l) . and so the observed similarity in the spec- 
tral type of the S-stars and HVSs is consistent with this scenario. Figure (HI shows an estimate of 



J HVS candidates are chosen for spectroscopy by color, to maximize the contrast against the halo popu- 
lation, and so are pre-selected to have B-type spectra (e.g. iBrown et a/...2006fii) . 
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the number of tidally captured S-stars for different MP populations (^Pe rets et'o'l ]l2007h . based on 
the observed orbital properties of young massive bina ries (a ^0. 2010 *^^ AU) and their fraction 
among young massive stars in the field (/2 ~ 0.75) dKobulnickv et al. l2006h . and on a model 
for the stellar density distribution in the in ner ^ 100 pc of the Galaxy (isothermal, normalized 
by the observations of Genzel et al. 2003h a nd a mass function model (continuous star forma- 
tion with a universal IMF, Figer et al. 12004 , see also Fig. |6]l. The typical binary was modeled 
as a 2 X 7.5 Af© binary (main sequence B-stars with a lifespan of ~ 5 x 10^ yr). Dynami- 
cal evaporation is negligible for such short-lived binaries. The steady state number of captured 
S-stars is then (A^*) = F/ci*. Figure (|4]i shows that with stellar relaxation alone, tidal capture can- 
not explain the S-star population. However, relaxation by GMCs is consistent with the observed 
number of S-stars, as well as with the spatial extent of the cluster of ~ 0.04 pc, which reflects 
the hardness of young massive binaries in the field (Eq. 12.2b . It is also consistent with the fact 
that the S-cluster does not include any star earlier then O8V/B0V. Such short-lived binaries are 
very rare in the field, and their mean number in the S-cluster is predicted to be {N^) < 1. 

The MP-induced binary tidal disruption scenario also predicts that there should be 10-50 
hyper-velocity ^ 4 Mq B-stars at distances betwee n 20 and 120 kpc from the GC. This is consis- 
tent with the total number of 43±31 extrapolated bv lBrown et al. (2006a, ). based on the HVSs de- 
tected at these distances in their field of search. The tidal disruption scenario predicts an isotropic 
distribution of HVSs around the GC, and a random ejection history, in co ntrast to mod els where 



the ejection is related to a discrete binary M BH merger event jYu & Tremaine 2003: Haardt et al. 
a iBaumgardt et al. 120061: iLevinI l2006h . The HVSs observed to date are consistent with an 



2006 



isotropic HVS distribution uniformly distributed in ejection time (iBrown et al.\\2006dt) and thus 
support the tidal disruption scenario. 

The tidal disruption scenario can naturally explain many of the properties of the S-stars and 
HVSs, but it has two potential flaws. (1) The predicted high eccentri cities of t he captured stars are 
larger than those observed for a few of the S-stars (e ^ 0.4. ,Eisenhauer et al 112005,) . However, the 
low observed eccentricities are expected to evolve after the capture by rapid resonant relaxation 
(® ■ (2) The lifespan of the most massive and shortest lived S-stars (i^, ~ 2 x 10*^ yr) is shorter 
by a factor < 10 than the MP-accelerated relaxation time in the inner ~ 5 pc (Fig. |2]), where a 
substantial fraction of the binaries are scattered from. Thus if a binary in those regions starts on 
a near-circular orbit, MP-induced relaxation is not fast enough to scatter it to a J < J/c orbit (Eq. 
11.2b within its lifetime. However, as the timescale discrepancy is not large, and as it affects only 
the most massive binaries in the central few pc, where the determination of Tr is ambiguous 
(see footnote on page|6]l, this does not appear to be a fatal flaw of this scenario. It does however 
highlight the importance of observationally quantifying the relaxation time in the GC and the 
distribution and properties of the field binaries. 

Low-mass bin aries are also deflec ted to the MBH by MPs and tidally disrupted at rates as high 
as ^ 10^^ yr~^ (iPerets et a/.ll2007l) . Neither the faint captured low-mass stars nor the late-type 
HVSs are detectable at this time. However, the captured stars affect the inner cusp dynamics in a 
way that may have implications for cosmic GW EMRI events. Binary disruption is effectively a 
local "source term" that modifies the flow of stars in phase space (cf Eq. [33] ). setting a diverging 
flow into the MBH and away from it, which modifies the steady state spatial distribution. Detailed 
calculations, which take into account the period distribution of low-mass binaries and the effects 
of binary evaporation, indicate that MP-induced tidal captures of white dwarfs close to the MBH 
efficiently competes against mass-segregation, which tends to lower the density of the low-mass 
white dwarfs there and raise the density of massive stellar BHs (§|3] Fig. [8]l. As a result, the 
cosmic rate of GW EMRI events involving white dwarfs is predicted to be at least comparable to 
that involving stellar BHs (Perets et al. 2007, in prep.). 

The proximity of the G C allows GW bursts from the fly-by of stars near the MBH to be de- 
tected (IRubbo et al ]l2006l) . MP-induced tidal binary disruptions increase the stellar density close 
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Figure 4. A comparison between the cumulative number of S-stars (main sequence B stars) observed or- 
biting the Galactic MBH on randomly oriented orbits (vertical bar), and the predicted number captured by 
3-body tidal interactions of the MBH with binarie s deflected to the center by massive perturbers, for differ- 
ent massive perturbers models l lPerets et a/.ll2007l) . The observed extent of the S-star cluster is indicated by 
the vertical hashed line. 



to the MBH and therefore the rate of GW bursts increases significantly. In partic ular, the rate of 
GW bursts from white dwarfs increases from ~ 0.1 yr~^ (iHopman et a/.ll2007h to a detectable 
rate of 2 yr^^ (Perets et al. 2007, in prep.). 

Binary MBHs form in the aftermath of galactic mergers, when the two MBHs sink by dy- 
namical friction to the center of the merged galaxy. Once the binary hardens, the orbital decay 
continues by 3-body interactions with stars that are deflected to the center and extract energy 
from the binary, until the orbit becomes tight enough for efficient GW emission, which rapidly 
leads to coalescence. Simulations show that when the loss-cone is replenished by stellar relax- 
ation alone^JhehUeractm rate is too slo w for the binary MBH to coales ce within a Hubble time 
(e.g. iBerczik et a/.ll2005l : see review by iMerritt & Milosavhevicll2005t Fig. |5]l. This "last par- 
sec stalling problem" appea rs to contra dict the circumstantial evidence that most galactic nuclei 
contain only a single MBH (IBerczik et a l. 2006; Merritt & Milosavljevic 2005| ), and furthermore 
implies few such very strong GW sources for LISA. On e rout for resolving the st alling problem 
is by accelerated MP-induced loss-cone replenishment (iPerets & Alexander l2007l) 

Figure (|5]l shows the time to coalescence, as function of the binary MBH mass, for different 
merger and MP scenarios, based on a combination of extrapolation of the Galactic MP population 
to early type galaxies, on extra-galactic observations of molecular gas in galactic centers, and on 
results from galactic merger simulations. The results show that MPs allow binary MBHs in gas- 
rich galaxies to coalesce within a Hubble time over nearly the entire range of M12. The situation 
with respect to gas-poor galaxies is less clear, since it is harder to model reliably the MPs there 



t Other p ossible routes are by interactions with gas in "w et mergers" (^anov ef cd. 1 19991: lEs cala et a/1 
20051: iDotti et al. 2QQ%, by interactions with a third MBH CMakino & Ebisuzakilll994l : iBlaes et al. 2002 
Iwasawa et al. 2006), or by accelerated loss-c one replenishment in a non-axisymmetric potential, ( 1YuI2003 : 
Berczik et a/.„2006) . or in a steep cusp llZieril2006h . 
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Figure 5. Accelerated binary MBH mergers in the presence of MPs iPerets & Alexan der 2007). Left: The 
time to coalescence as function of binary MBH mass, for different merger scenarios distinguished by the 
mass ratio Q between the two MBHs and the MP contents of host galaxies. The age of the universe is 
indicated by the dotted horizontal line. Stellar relaxation alone cannot supply a high enough rate of stars 
for the slingshot mechanism to complete the merger within a Hubble time. However, in minor mergers 
(Q = 0.05) major gas-rich mergers {Q = 1) with MPs merger is possible within a Hubble time for all but 
the most massive MBHs. Right: The evolution of the binary MBH semi-major axis as function of time for 
major mergers (Q — 1) in the presence of MPs (solid line) and stellar relaxation alone (dashed line), for 
binary MBH masses of 10®, 10^, 10** and lO'' Mq (from bottom up). 



(probably clusters rather than GMCs). However, even for such galaxies, MPs allow coalescence 
within a Hubble time up to masses of A/12 10* ^^o- 

Efficient binary MBH coalescence by MPs has various implications. It increases the cos- 
mic rate of GW events from MBH-MBH mergers, it increase the "mass d eficit" in the galac 



tic core (the steUar mass ejected from the core by the slingshot ef fect) ( IMilosavlievic et al. 



2OO2I: iRavindranath et a/.ll2002t lOrahaml I2OO4I: iFerrarese et a/.ll2006l) . it leads to the ejection 



of hyper-velocity stars to the inter-galactic space, but it s uppress the formation of triple MBH 
systems and the ejection of MBHs into interg alactic space ( Saslaw et aZ.ll974 : Blaes et aZll2002 : 
iHoffman & LoeblboOTt llwasawa et a/.ll2006h . 



3. Strong mass segregation 

3.1. The Bahcall-Wolf solution of moderate mass-segregation 

The 2-body relaxation timescale around the Galactic MBH, Tr ^ 0(1 Gyr), is short enough 
for the old stellar population there to relax to a universal steady-state configuration, indepen- 
dently of the initial conditions. This configuration was investigated by Bahcall & Wolf (197(^ 



1977h . The Bahcall-Wolf solution predicts that in the Keplerian potential near a MBH, stars 



of mass in a multi-mass population ,Mi < < M2, have a DF that is approximately a 
power-law of the specific orbital energy e, /m (e) oc t^" , where pm oc A/^ with a proportionality 
constant Pm/M^, ~ l/(4Af2)- In a Keplerian potential, this DF corresponds to a density cusp 
JT-Af ocr""*^ , where a a/ = 3/2 + pm- Elementary consid erations show that q: = 7/4 (p= 1/4) 



for a single mass population (e.g. iBinnev & Tremainell987 , § 8.4-7). This follows from the con 



servation of the orbital energy that is extracted from stars that are scattered into the MBH, and 
transferred outward by the ambient scattering stars in a steady-state, distance-independent cur- 
rent, dE{r) /dt ^Ni,{<r)Ei, (r) /Tr (r) ~ r'^/^"^" = const (using the relations A'^* (< r) cx r^"" , 
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Figure 6. Left: A theoretical population model for the central few pc of the GC llAlexander & Sternberg 
[T999; Alexander 200 ^, assuming cont inuous star formation over the past 10 Gyr jFiger et fl/.ll2004l) with a 
"universal" IMF ( Miller & Scalolll979h . Bottom panel: The fraction of old stars (defined here as stars with 
main-sequence lifespa n of > 1 Gyr) in the pop ulation, as function of the iv -band magnitude in the GC (for 
DM+Ak = 17.2 mag. lEisenhauer et fl/.l200a) . Middle panel: The JiT-band luminosity function. Top panel: 
The mean mass of all stars and of the old stars only, as function of the A'-band magnitude. The concentration 
of old Red Clump / horizontal branch giants around ~ 15.5 is clearly seen as an excess in the luminos- 
ity function, as an increase the fraction of o ld stars and as a dec rease in the mean stellar mass relative to 
stars both immediately brighter and fainter fechodel et alJ^lOO'f). Right: The observed azimuthally-aver- 
aged stellar surface number density around the Galactic MBH as function of projected angular distance, 
E(7i) (R = 1" corresponds to ~ 0.04 pc in the GC) in 3 adjacent J\ -magnitude bins, centered around the 
bin associated with the Red Clump giants (14.75 < K < 15.75), with broken power-law fits E cx 
dSchodel et al\2QffK . adapted with permission from Astronomy and Astrophysics). Top: 15.75 < AT < 16.75 
(density multiplied by 2 for display purposes). Middle: 14.75 < K < 15.75 (the Red Clump / horizontal 
branch range). Bottom: 13.75 <K < 14.75. 



Ei, ~ and Tr oc ^"-3/2^ ^1.21 ). The Bahcall-Wolf solution reproduces this result for a single 
mass population, and predicts that it should apply also to the heaviest stars in a multi-mass pop- 
ulation. The Bahcall-Wolf solution thus implies that at most Aa =1/4 between the lightest and 
heaviest stars in the population. The predicted degree of segregation is moderate. 

Theoretical considerations, results from dynamical simulations and GC observations, hint that 
the moderate segregation solution should not and does not always hold, even in relaxed systems. 
As formulated, the solution depends only on the stellar masses, but not on the mass function. 
However, this cannot apply generally, since in the limit where the massive objects are very rare, 
they are expected to sink efficiently to the center by dynamical friction, and create a cusp much 
steeper than a = 7/4. As shown below 03.2I ). models of the present-day mass function in the cen- 
tral few pc of the GC suggest that the massive objects are relatively rare. Dynamical simulations 
of mass segregation in the GC based on such a mass function (Fig. [8]) indeed show steep cusps 
(a > 2) for the heaviest masses. Finally, the observed surface density distribution of GC stars in 
the magnitude bin 14.75 < K < 15.75, which corresponds to the low-mass (0.5 < M^, < 2 Mq) 
Red Clump / horizontal branch giants (Fig.|6|, is substantially flatter than that of the higher-mass 
giants (M^ ^ 3 Mq ) that populate the adjacent bins of brighter and fainter magnitudes (Fig.|6l 
Schodel et al. .l2007h . The sign of this trend is as expected for mass segregation, but the size of the 



effect is much larger than predicted by the Bahcall-Wolf mod erate segreg ation solution. How- 
ever, it can be explained in terms of mass-segregation if Aa > 1 ( Levil2006 ). While none of these 
hints for strong mass-segregation is decisive in itself, and other explanations are possible, taken 
together they motivate a re-examination of the mass-segregation solution in a relaxed system. 
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3.2. The relaxational self-coupling parameter 

Assume for simplicity a stellar system with a two-mass population of light stars of mass Af^, 
total initial number Nl and local density n l (r) and heavy stars of mass AIh, total initial number 
Nh and local density nn (r) . The self interaction rate is then Fli, oc ulM^ /v^ for the light stars 
and Thh oc uhM'^/v^ for the heavy stars (©. In the limit where the heavy stars are test 
particles {nn/nL ^ Mf^/Mfj, or equivalently Thh /^ll ^ 1), the heavy stars interact mostly 
with the light ones, lose energy and sink to the center by dynamical friction. Conversely, in the 
limit Thh IT ll ^ 1, the heavy stars interact mostly with each other, effectively decouple from 
the light stars and establish an a = 7/4 cusp typical of a single mass population. This suggests 
that the global relaxational self-coupling parameter (cf Eq. 12.1b . defined as 



P2 = NhMI/NlMI, (3.1) 

can be used to determine whether the system settles into the moderate (Bahcall-Wolf) mass- 
segregation solution (/i2 > 1) or the strong mass-segregation solution (/i2 < 1). This hypothesis 
is borne out by the numerical results presented below|j (Alexander & Hopman 2007, in prep.; 
93.3b . For a continuous mass distribution, ^2 can be generalized to 

P2= M^{dN/dAh)dM^ / / M^idN/dAQdM^, (3.2) 

the ratio between the 2nd moments of the mass distribution of the heavy [M^, > Mq) and light 
{Mi^ < Mq) stars, for some suitable choice of the light/heavy boundary mass Mq. 

The value of p2 depends on the population's present-day mass function. So-called universal 
initial mass functions (IMFs), which extend a l l the w ay from the brown dwarf boundary Mi ~ 



0.1 Mp, to M2 > 100 M fTj (e.g. the ISalpeted (Il955h IMF, and its subsequent refinements, the 
Miller & ScalolJl979l) and lKroupal (1200 ll) IMFs), result in evolved populations, old star-bursts or 
continuously star forming populations, that naturally separate into two mass scales, the 0(1 Mq) 
scale of low-mass main-sequence dwarfs, white dwarfs and neutrons stars, and the 0(10 Mq) 
scale of stellar black holes, and typically have < 1 (Fig.Q. Such evolved populations are thus 
well-approximated by the simple 2-mass population model. In particular, the volume-averaged 
stellar population in the central few pc of the GC is reasonably well approximated by a 10 Gyr 
old, continuously star-forming population with a universal IMlQ dAlexander & Sternberg|[l999 



Fig. |6]l. Generally, 10 Gyr old, continuously star-forming populations with a power-law IMF, 
dN/dM-i, cx M^^, have /i2 < 1 for 7 > 2, and fi2> 1 for 7^2. Since the critical value 7 = 2 is 
close to the generic Salpeter index 7 = 2.35, it is quite possible that both the moderate and strong 
segregation solutions are realized around galactic MBHs, depending on the system-to-system 
scatter in the IMF (and perhaps also realized in clusters around IMBHs, if such exist). 



3.3. Solutions of the Fokker-Planck energy equation 

The steady state configuration of stars around a MBH can be described in terms of the diffusion 
of stars in phase space, from an infinite reservoir of unbound stars with a given mass function 
(the host galaxy, far from the MBH), to an absorbing boundary at high energy where stars are 
destroyed (the MBH event horizon, tidal disruption radius, or collisional destruction radius). 

f In the limit Mh /M l 2> 1, it may be necessary to take explicitly into account the the dynamical friction 
timescale in order to obtain a more accurate segregation criterion. Here fj,2 is adopted for its simplicity. 

X Note that recent analysis of late type giants in the GC suggests tha t the IMF in the inner ~ 1 pc of 
the GC could typically be a flat 7 ~ 0.85 power-law jManess et fl/.ll2007l) . This would imply ^ 1 in the 
inner ~ 1 p c, possibly a volume-averaged > 1 in the inner few pc (the "collection basin" for stellar BHs, 
iMiralda-^cude & Gould,2000.) . and hence moderate segregation. 
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Figure 7. The predi cted values of the global relaxational self-coupling parameter ^2 for a "universal" 
iMiller & Scald l ll979h IMF (top line, shifted by x 100 for display purposes), an evolved mass function 
assuming continuous stars formation over 10 Gyr (middle line), and an evolved star formation burst 10 Gyr 
old (bottom line, shifted by xO.Ol for display purposes) (Alexander & Hopman 2007, in prep.). The mass 
functions of the old populations develop excesses in the ~ 0.6-1.4 Mq range due to the accumulation of 
white dwarfs and neutron stars, and in the ~ 10 Mq range due to the accum ulation of stellar black holes 
(here represented by a simplified discrete mass spectrum, see I Alexande j|200^ table 2.1). 



Bahcall & Wol(ll976ifl977l) simplified the full Fokker- Planck treatment in (E, J) phase space 



by integrating over J so as to reduce it to E only, by assuming a Keplerian potential, and by 
recasting i t in the form of a particle con servation equation. In dimensionless form these can be 
written as dHopman & Alexandeill2006fcl) 



d d 

■7r9Mix, t) = -x^^^—QAiix, t) - Rm{x, t) , (3.3) 
OT ox 

where M, x and r are the dimensionless mass, energy, and time, respectively, gj\/ is the di- 
mensionless DF, Q M is the flow integral, which expresses the diffusion rate of stars by 2-body 
scattering to energies above x, and Rm oc /Tr is the J-averaged effective loss-cone term. Qm 
and are non-linear functions of the set of DFs {.gj\/}. The equations are solved for by 
finite difference methods starting from an arbitrary initial DF and integrating forward in time un- 
til steady state is reached, subject to the boundary conditions that no stars exist at energies above 
some destruction energy xjj, gM{x > xd) = 0, and that the unbound stars ar e drawn from an 
isother mal distribution with a given mass function, (7a/(x' < 0) = cxp(Afx). lBahcall & WoH 



(11977b showed that the stellar space density distribution. 



nM{r)oc / gM{x)^r/rh ~ xAx , (3.4) 



does not depend strongly on the exact form of the loss-cone term, and proceeded to use in their 
mass-segregation calculations a simplified version of Eq. ( I3.3l l by setting Rm ~ 0. This approx- 
imation can be justified by noting that while the existence of a loss-cone drastically increases 
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Figure 8. Numerical models of the mass distribution in the GC showing strong segregation. Right: The 
density distribution of 10 A/© stellar BHs and 0.7 A/q white dwarfs at 10 Gyr in an approximate model of 
the GC, with the evol ved universal IMF o f Fig. derived by M. Freitag (priv. comm., reproduced here with 
permission. See also iFreitag et al using an implementation of the Henon method (Freitag & BenJ 

[2002). The low -mass white dwarfs settle into a ol ~ 1.5 power-law cusp. The distribution of massive 
stellar BHs can be approximated by a piece-wise broken power-law, with oih ~ 2.2 at r ~ 0.1 pc. Left: 
The space density of a simplified 4-component population model fo r the GC, as given by the so lution 
of the Fokker-Planck equations with a loss-cone term (adapted from lHopman & Alexandeill2006/i with 
permission from the Astrophysical Journal). The logarithmic slope of the density cusp of the stellar BHs at 
r = 0.1 pc is a~ 2.1, as compared to a ~ 1.5 for the lighter species. 



the flow rate of stars into the MBH, it typicafly affects only a small volume in phase space 
near J ^ Q. This translates to small changes only in the J-integrated DF gM{x), mostly for 
X xd, and even smaller changes in njv/(r) due to the smoothing effect of the gM{x)^nM(,r) 
transformation (Eg. 13.4b . Here we adopt this approximation to allow direct comparison with the 



Bahcall & Wolf! ( Il977h results, after verifying that the inclusion of the loss-cone term indeed 



does not significantly change the derived stellar cusps (cf Fig.[8]and|9]l. 

We calculated a suite of such Fokker-Planck mass-segregation models for 2-mass popula- 
tions with different mass ratios Mh/Ml and mass functions Nh/Nl, spanning a very wide 
range of the global relaxational self-coupling parameter valueqj, 10^"^ < M2 < 10'^ (Alexan- 
der & Hopman 2007, in prep.). The DFs are not exact power-laws, and the logarithmic slopes 
Pm {x) = d loggA/ / d logx depend somewhat on energy, especially near the boundaries. However, 
analysis of the results is considerably simplified by the fact that the values of pai{x) vary mono- 
tonically with /i2 at all x, and so the order ranking of pM for different models does not depend 
on the choice of x. Figure (|9|l shows and pn at x ~ 10, which corresponds to r ~' 0.1 pc 
in the GC. This choice samples gAi{x) in a representative region, far from either boundaries at 
a; = and xd ~ 10^, and translates to an observationally relevant region in the GC, which is close 
enough to the MBH to be nearly Keplerian, but still contains a large number of observed stars to 
allow meaningful statistics (cf Fig.|6]l. 

Figure ^ shows that for /i2 > 1, the Fokker-Planck calculations recover the Bahcall-Wolf 
solution: pn ~ 1/4 irrespective of the mass ratio, andp^ ~ {1/4){Ml/AIh)- However, for /.(2 < 1 
there is a marked qualitative change in the nature of the solutions, as anticipated by the analysis in 
33.21 The more the light stars dominate the population (the smaller fi2), the more they approach 
the single population solution = 1/4 ( 33.11 ). The heavy stars s ettle to a much steeper c usp with 
Ph > 1/4. Figure (|9]l also shows the grid of models explored by Bahcall & Wolf! ( 1977 ). which, 
while large, covers only the /i2 > 1 range. This explains why the strong segregation branch of the 



f It is unlikely that real stellar system will have relaxational self-coupling parameters /i2 <C 0.1. How- 
ever, the study of such models is useful for understanding the mathematical properties of the solution. 
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solutions escaped their notice (their one model with /^J2 ~0.6 has a low mass ratio Mh/AIl = 1.5, 
where the two solution branches are not very different). 

The ^2 < 1 models explored here follow the pj^j oc Ali, relation noted by iBahcall & WoH 
(Il977h for the approximate mass-segregation solutions without a loss-cone term. Therefore, in 
those models where the limit ^ 1/4 is reached (for Mh/Ml = 1.5, 3), the heavy stars reach 
the asymptotic value p//^ {\/A){Mh /Ml)- It remains to be seen whether this result also holds 
for higher mass ratios, and for the full Fokker- Planck equation (Eg. 13.3b with the loss-cone term. 

A realistic evolved stellar system, such as the GC, is expected to have a maximal mass ratio 
of at least Mh/Ml = 10 and ^2 — 0.15 (Fig. |7]i. The mass-segregation calculations indicate 
that for these parameters the stellar BHs are expected to form an an — 2.1-2.2 cusp (Figs. [8] 
|9]l, significantly steeper than the an = 1.75 predicted by the Bahcall-Wolf solution of moderate 
segregation. It is encouraging that this logarithmic slope is very close to that found in numerical 
simulations (Fig. |8]l and that it is broadly consistent with what is needed to explain the observed 
trend in the stellar surface density distributions in the GC in terms of mass segregation (Fig.|6]l. In 
other systems the moderate s egregation solution m ay apply. For example, if the globular cluster 
Ml 5 contains an I MBH (e.g. iGerssen et al 1 l2002l) . then a tentative determination of its present- 
day mass function ( Murphv ef a/.lll997l) suggests a high relaxational self-coupling parameter, 
/X2 ~ 40 and a relatively shallow a = 7/4 cusp of stellar BHs. Full-scale numeric simulations that 
are free of the restrictive assumptions of the analytic approach adopted here (Keplerian potential, 
fixed boundary conditions, approximate treatment of the loss-cone and a fixed 2-mass stellar 
population) are needed to verify and test these predictions in more detail. 

Strong segregation will affect the cosmic rates of EMRI events. A detailed analysis of the 
anticipated change relative to the various discrepant published rate estimates depends on their 
specific assumptions (e.g. the assumed mass function, slope of the cusp, normalization of the 
stellar number density), and is outside the scope of this review. 



4. Resonant relaxation 

4. 1 . Resonant relaxation dynamics 

The effect of 2-body relaxation on a test star is incoherent: the star experiences randomly ori- 
ented, uncorrected perturbations from the ambient stars, and as a result its orbit deviates in a 
random-walk fashion from its original phase-space coordinates (in a stationary spherical smoothed 
potential where E and J would have been conserved in the contin uum limit, A_E oc and 
A J C)C \/t due to 2-b ody interactions). Resonant relaxation (RR) (IRauch & Tremaine 119961: 



Rauch & Ingallslll998h is a form of accelerated relaxation of the orbital angular momentum, 
which occurs when approximate symmetries in the potential restrict the orbital evolution of the 
perturbing stars. This happens in the almost Keplerian potential near a MBH, where the orbits 
are approximately fixed ellipses (the potential of the enclosed stellar mass far from the MBH, or 
General Relativistic (GR) precession near the MBH, eventually leads to deviations from pure el- 
liptical orbits), or in a non-Keplerian, but nearly spherically symmetric potential around a MBH, 
where the orbits approximately conserve their angular momentum and move on rosette-like pla- 
nar orbits (the fluctuations in the potential due to stellar motions eventually lead to deviations 
from strictly planar orbits). As long as the symmetry is approximately conserved, on times shorter 
than the coherence timescale t^j, the orbit of a test star with semi-major axis a experiences cor- 
related (coherent) perturbation^ which can be described as a constant residual torque exerted 
by the superposed potentials of the N-i,{< a) randomly oriented elliptical "mass wires" (in a 
Keplerian potential ) or "mass annuli" (in a non-Keplerian spherical potential) that represent the 



f RR is better described as "coherent relaxation". The term "resonant" refers to the equality of the radial 
and azimuthal orbital periods in a Keplerian potential, which results in closed ellipse orbits. 
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Figure 9. Fokker-Planck mass-segregation results. The logarithmic slopes pu and pl of the DFs of the 
heavy stars (thick hnes) and light stars (narrow lines) , evaluated at (r ~ 0.1 pc in the GC), as function of 
the global relaxational self-coupling parameter /i2, for mass ratios of Mh /Ml = 1.5, 3, 10 (Alexander & 
Hopman 2007, in prep.). The logarithmic slope of the stellar density of massive stars in the GC, assuming a 
universal IMF and continuous star formation history (^2 ~ 0.13, Fig.|7] = 3/2 + ph — 2.2) is indicated 
by a cross on the left, and for globular cluster Ml 5 ( assuming it harbors an IMBH), on the right(estimated 
at ^2 ~ 37, based o n the mass function mo del of iMurphv et al\ ( 119971) , ajj — 1.75). The results for the 
models calculated bv lBahcall & Wol3 ( 1 19771) are indicated by circles. 



orbitally-averaged mass distribution of individual perturbing stars. The magnitude of the residual 

1 /2 

torque is then J ~ iV^ (< a)G'A/*/a and the change in the angular momentum of the test star 
increases linearly with time, AJ^Jt (for t < t^). The orbital energy, on the other hand, remains 
unchanged, since the potential is constant. 

RR in a Keplerian potential is called scalar RR since it changes both the magnitude and direc- 
tion of J. Scalar RR can therefore change a circular orbit into an almost radial, MBH-approaching 
one. In contrast, RR in a non-Keplerian spherical potential is called vector RR since, for reasons 
of symmetry, it changes only the direction of J, but not its magnitude (Fig.fTOli. Vector RR can 
randomize the orbital orientations, but does not play a role in supplying stars to the loss-cone. 

On timescales longer than the coherence time, the orbital orientations of the perturbing stars 
drift, and coherence is lost, the maximal change in angular momentum during the linear coher- 
ence time, AJ^ ^ Jt^ then becomes the "mean free path" for a random walk in J-space, whose 
time-step is t^. On timescales longer than the coherence time, the angular momentum changes 
incoherently cx ^/t, but much faster than it would have in the absence of RR. The energy is 
unaffected by RR and always evolves incoherently cx on the long non-resonant relaxation 
timescale (Fig. fTTT i. The RR timescale Trr is defined, like the incoherent 2-body relaxation 
timescale, as the time to change J by order Jc (Eqs. ll.llll.2b . Trr^ (Jc/AJ^j)^ C which can 
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Stationary ellipses 
in point mass potential 



Planar rosettes in 
spherical potential 





Scalar resonant relaxation 





Vector resonant relaxation 



Figure 10. A sketch comparing the symmetries leading to scalar and vector RR. Top: The torques by fixed 
elliptical "mass wires" in a Keplerian potential lead to rapid changes in both the direction and magnitude of 
the orbital angular momentum of a test star. Bottom: The torques by fixed "mass annuli" in a non-Keplerian 
spherical potential lead to rapid changes in the direction, but not in the magnitude of the orbital angular 
momentum of a test star. 



be expressed as ( Hopman & Alexandedl2006al) 

N.{<a) P\a) 



Trm = A 



A 



RR 



RR 



^2(< a) U 



N4< a) 



Ah 



(4.1) 

N^A'Q, P is the 



where fj. is the relative enclosed enclosed stellar mass, fj, ~ N^Mi,/ (A/, 
radial orbital period, and the last approximate equality holds in the Keplerian regime. Here and 
below, the constants A*^^ are numerical factors of order unity that depend on the specifics of the 
coherence-limiting process, on the orbital characteristics of the test star, and probably also on the 
parameters of the stellar distribution. These constants are not well-determined at this time. 

The coherence time depends on the symmetry assumed and on the process that breaks it. For 
a non-relativistic near-Keplerian potential, the Umiting process is precession due to the potential 
of the distributed stellar mas^ 



M, 



P{a). 



(4.2) 



Remarkably, the resulting RR timescale does not depend on N^,. Close to the MBH is much 
shorter than the non-coherent 2-body relaxation timescale (here denoted for emphasis as T^r), 



T. 



M 
RR 



N4< a)Ah 



T, 



NR. ■ 



(4.3) 



Yet closer to the MBH, it is GR precession that limits the coherence. 



toR 



8 / J 

3 V'^LSO 



Pia). 



(4.4) 



where Jlso = 4:GM,/c is the last stable orbit for e ^ c^. The GR precession is prograde, while 

f The enclosed stellar mass iV*A/* changes the Keplerian period PocA/. by AP/P^ N^Ah/2M, 
= Aip/2TT. Identifying de-coherence with a phase drift Aip — TT then implies f^/ ~ {n/ Aip)P. 
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Figure 1 1 . A correlation analysis of A'-body simulations (A = 200) showing the relaxation of energy and 
of scalar and vector angular momentum around a MBH in the Keplerian limit (Af*/M, = 3 x 10~^) for a 
thermal population of stars, as function of the elapsed time-lag (Eilon, Kupi & Alexander, 2007, in prep.). 
The change in /S.E/E, AJ/ Jc and | AJ| / Jc is plotted as function of the normalized time-lag r = At/P 
in the range 1 < r < 100. The mass precession coherence time of the system is tm = tM/ P 1-7 x 10^, 
and the potential fluctuation coherence time is T^ = t^/ P ^ 1.2 x 10'', so both scalar and vector RR are 
expected to grow linearly over the plotted range. More detailed analysis shows that the AJ/ Jc(E) is a 
function of both energy and angular momentum, which for t ^ scales as ^/r, and for 1<S^t<^Tw scales 
as T, and that | AJ| / Jc is simply proportiona l to A J/ Jc The correlatio n analysis is an efficient method for 
quantifying relaxation in A-body results (cf lRauch & Tremainelll99^ , Fig. 1). The theoretical predictions 
for A J/ Jc and | AJ| / Jc fit the numeric results very well. As expected, AE / E(x^/t at all time-lags. 

that due to the distributed mass is retrograde, and so they may partially cancel each other Their 
combined effect on the scalar RR timescale is 

2 



A 



RR 



RR 



N,i< a) 



1 



1 

^GR 



(4.5) 



Since tM increases with r, while tcR decreases with r, scalar RR is fastest at some finite distance 
from the MBH, which typically coincides with ~rcrit/2 for LISA EMRl targets (Fig.[T3]l. 

Precession does not affect vector RR. The coherence in a non-Keplerian spherical potential is 
limited by the change in the total gravitational potential (f>~4>, + (j)^, caused by the fluctuations 
in the stellar potential, 0^, due to the realignment of the stars as they rotate by tt on their orbits. 
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where the last approximate equality holds in the Keplerian regime. The vector RR timescale is 
obtained by substituting in Eq. ( 14.1b . 



rpV 

^RR 
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Vector RR is much faster than scalar RR (Fig.fTSll. 
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4.2. Resonant relaxation and EMRI rates 

The efficiency of scalar RR quickly decreases with distance from the MBH, since the coher- 
ence time falls as Af^(< r)/M, grows. At rh, where M^(< ry,) /M, ~ 0(1), scalar RR is 



almost completely quenched. Because Tcrit ~ fh for tidal disruptio n dLightman & S hapiro 1 977 



91. U . RR does not significantly enhance the tidal disruption rate dRauch & Tremaineiil996 ). In 



contrast, rcrit ^ 0.01 pc for EMRI events, where M,^(< rcrit)/Af« ^ 1 and T^^ is near its 
minimum. Scalar RR therefor e dominates the dynamics of the loss-cone for GW EMRI events 
( Hopman & Alexanderl2006a ). Scalar RR accelerates the flow of stars in phase-space from large- 



J orbits to low- J orbits that approach the MBH and can lose orbital energy and angular momen- 
tum by the emission of GWs. However, if unchecked, RR would continue to rapidly drive the 
stars to plunging orbits that fal l directly into the MBH. This is where GR precession is pre- 
dicted to play an important role ( Hopman & Alexandej2006a ). Orbits with very small periapse. 



Tp ~ fewXTs, where GW emission becomes appreciable, are also orbits where the GR precession 
rate becomes large enough (tcR becomes short enough) to quench RR, and allow the EMRI in- 
spiral to proceed undisturbed. This subtle cancellation, which is critical for the observability of 
EMRI events, still has to be verified by direct simulations. 

The effect of scalar RR can be included in an approximate way in the Fokker-Planck equa- 
tion (Eq. I3.3l l as an additional loss-cone term Rrb. oc Xd/T^p, where the efficiency factor x 
parametrizes the uncertainties that enter through the various order-unity factors Af^^j^ (Eq. |4.1| l. 
Such calculations show that the poorly determined value of the efficiency can strongly affect the 
predicted EMRI rates (Fig. [T2b . As the efficiency rises, the EMRI rate first increases because 
stars are supplied faster to the loss-cone, but when the efficiency continues to rise the stars are 
dr ained so rapidly into the M BH, that the EMRI rates are strongly suppressed. 



Rauch & Tremaine (Il996h explored the efficiency of RR by a few small-scale iV-body simu- 



lations, and noted a large variance around the derived mean efficiency. Here we use their mean 
efficiency as the reference point (x = 1), but consider also values smaller and larger. Figure 
ST% shows the GW inspiral rate and direct plunge rate as function of the unknown efficiency 
X, relative to the no-RR case (x = 0), derived from Fokker-Planck calculations of a single mass 
population. The EMRI rate rises to ^ 8 times more than is expected without RR, peaking at 
1 < X < 2, but then falls rapidly to zero at x^t 10- The strong x-dependence of the EMRI rates 
provides strong motivation to deter mine the RR efficiency an d its dependence on the parame- 



ters of the system both numerically jGiirkan & Hopmanll2007l : Eilon, Kupi & Alexander 2007 



in prep; Fig [TTI ). and by direct observations of the only accessible system at present where RR 
effects may play a role — the stars around the Galactic MBH. 

4.3. Resonant relaxation and stellar populations in the GC 

The stellar population in the GC includes both young and old stars, and is composed of distinct 
sub-populations, each with its own kinematical properties (see Alexander 200i for a review). As 
shown below, RR can naturally explain some of the systematic differences between the various 
dynamical components in the GC. Conversely, GC observations of these populations can then 
test the various assumptions and approximations that enter into analytic treatment of RR, and in 
particular constrain the poorly determined RR efficiency. 

Figure ( fT3] ) summarizes the typical distance scales and ages associated with these populations, 
and compares them with the various relaxation timescales. The calculation of the relaxation 
timescales are approximate since they assume a single mass population. The non-resonant 2- 
body relaxation timescale (Eq. ll.lb is roughly independent of radius in the GC, T/vij ^ few x 10^ 
yr (assuming a mean stellar mass of A/* = 1 M©; in a multi-mass s ystem it is expected to decreas e 



to T^vfl ~ 10 yr in the inner 0.001 pc due to mass segregation. iHopman & Alexandeill2006d) . 
Because neither the RR efficiency, nor the mass function is known with confidence, the scalar 
RR timescale, T^j^, is shown for two different assumptions; X-^^* = 1 ^& and X-^* = 10 Mq. 




Rauch & Tremaine 1996 



10 ^ 10 ' lO" lo' 

Resonant relaxation efficiency 

Figure 12. The relative rates of GW EMRl events and direct infall (plunge) events, as function of the 
unknown efficiency of RR, x, normalized to x = 1 for the values derived by .Rauch & Tremaine ( 199^ . 



As discussed in ^4.11 beyond r 0.1 pc, T|,^ > Tjsir due to mass precession, and the loss- 
cone replenishment is dominated by non-coherent relaxation. T|,^ decreases toward the MBH, 
until it reaches a minimum, where it starts increasing again due to GR precession. The distance 
scale where T^^ is shortest happens to coincide with the volume r < r„it, where most GW 
EMRI sources are expected to lie and where <C T/v i?> so RR dominates EMRI loss-cone 
dynamics 04.2b . In contrast to scalar RR, the vector RR timescale (shown here for an 
assumed x^^* = 1 -^0) decreases unquenched toward the MBH. 

Dynamical populations and structures whose estimated age exceed these relaxation timescales 
must be relaxed. Those whose age cannot be determined, but whose lifespan exceeds the relax- 
ation timescales may be affected, unless we are observing them at an atypical time soon after 
they were created. The younges t dynamical structure observed in the G C is the stella r disk (or 
possib ly two non-aligned disks) ( Levin & Beloborodovl2003 : Genzel et a l. 2003; Pau mard et al. 
20061) . which is composed of ~ 50 young massive OB stars with an age of t^, ~ 6 ± 2 Myr, on 



co-planar, co-rotating orbits that extend between ^ 0.04-0.5 pc. The inner edge of the disk is 
sharply defined and it coincides with the outer boundary of the S-stars cluster 02.2b . Figure ( flj]) 
shows that the vector RR timescale equals the age of the stellar disk at its inner edge, and so is 
consistent with the spatial extent of the disk. Even if the S-stars were initially the inner part of 
the disk (this does not appear likely given that they are systematically lighter than the disk stars), 
vector RR would h ave efficiently randomized their orbital inc lination. However, their measured 
high eccentricities ( lEisenhauer et a/.ll2005l: iGhez et a/.ll2005b would then be hard to explain. If 
instead, the S-stars are not-related to the disks, but were tidally captured around the MBH by 
3-body exchange interactions 02.31 1, then only their lifespan can be determined. Tidal capture 
leads to an extremely eccentric captured orbit (Eq. l2.2l i. Scalar RR could then randomize and de- 
crease the eccentricities of at least a few of the older S-stars closer to the MBH. Vector and scalar 
RR could also exp lain why the old evolved giants (with pro genitor masses of M j, ^ 2-8 Mq, 
Genzel et al. I I1994 at r > 0.1 pc appear dynamically relaxed (iGenzel et al. boooh . in spite of the 



fact that their lifespans are shorter than the non-coherent relaxation time. 
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Figure 13. Evidence for resonant relaxation in the GC in the age .vs. distance from the MBH plane. 
The spatial extent and estimated age of the various dynamical sub-populations in the GC (shaded ar- 
eas) is compared with the non-resonant 2-body relaxation timescale (top line, for assumed mean mass 
of Af* = 1 Mq) and with the scalar RR timescale (two curved lines, top one for x^^* = 1 Mq, bottom one 
for x^* = 10 Mq) and vector RR timescale (bottom line, for xM* = 1 Mq). The populations include the 
young stellar rings in the GC (filled rectangle in the bottom right); the S-stars, if they were bom with the 
disks (open rectangle in the bottom left); the maximal lifespan of the S-stars (filled rectangle in the middle 
left); the dynamically relaxed old red giants (filled rectangle in the top right); and the reservoir of GW 
inspiral sources, where the age is roughly estimated by the progenitor's age or the time to sink to the center 
(open rectan gle in the top left). Stellar com ponents that are older than the various relaxation times must be 
randomized. jHopman & Alexandeill20063 . reproduced with permission from the Astrophysical Journal). 

It should be noted that the effect of RR on the stellar density distribution is not expected to be 
large even quite close to the MBH (r < 0.1 pc), unless the efficiency x is very high, because the 
RR-induced depletion of the DF at high energies is smoothed by the transformation from the DF 
to n*(r) (Eq. l3.4l i and by the contribution of unbound stars to the central density. 

5. Summary 

Relaxation processes play an important role in the GC, where the 2-body relaxation time is 
shorter than the age of the system and the stellar density is high. The scaling laws that follow 
from the M, / a relation imply that the same must hold for all galaxies with A/, < few x 1 0^ Mq . 
Relaxation processes affect the distribution of stars and compact remnants, lead to close interac- 
tions between them and the MBH, and may be related to the unusual stellar populations that are 
observed in the GC. These are of relevance because of the very high quality stellar data coming 
from the GC, and because galactic nuclei with low-mass MBHs like the GC are expected to be 
important GW EMRl targets for the next generation of space borne GW detectors. In addition, 
efficient relaxation mechanisms that operate and can be studied in the GC may play a role even 
in galactic nuclei with high-mass MBHs, where 2-body relaxation is unimportant. 

Three processes beyond minimal two-body relaxation were discussed here: accelerated loss- 
cone replenishment by MPs, strong mass-segregation in evolved populations, and rapid RR. Ev- 
idence was presented that these processes operate and may even dominate relaxation and its 
consequences in the GC: The S-stars and HVSs are consistent with relaxation by GMCs; there 
are hints for strong mass segregation in the central density suppression of the low-mass Red 
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Clump giants and in numeric simulations of the GC, and RR appears to play a role in the trun- 
cation of the stellar disks and the orbital randomization of the S-stars and the late type giants. 
There are also cosmic implications: MPs enable the efficient merger of binary MBHs, and boost 
the rates of white dwarf EMRIs captured near the MBH by tidal disruptions of stellar binaries. 
Strong segregation, and in particular RR can strongly affect the EMRI rates from stellar BHs. 

The stellar dynamics laboratory in the GC holds great promise for future progress in under- 
standing these mechanisms and their implications. 
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